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This paper presents a general investigation into the improvement of modal scaling 
factors of an experimental modal model using mass additive technique. Data base 
required by the proposed method consists of an experimental modal model (a set of 
complex eigenvalues and eigenvectors) of the original structure and a corresponding 
set of complex eigenvalues of the mass-added structure. Three analytical methods, 

i.e., first order and second order perturbation methods, and local eigenvalue 
modification technique, are proposed to predict the improved modal scaling factors. 
Difficulties encountered in scaling closely spaced modes are discussed. Methods to 
compute the necessary rotational modal vectors at the mass additive points are also 
proposed to increase the accuracy of the analytical prediction. 


INTRODUCTION 


For most applications of experimental modal database, the identified modal vectors are expected to be 
normalized according to {y r } r [ U K^r} = 1 • But in practice, this relation is rarely satisfied and it becomes 
O' r } T [ U]{y r } = Va? , where a T is an unknown scaler of the r 01 modal vector. 

Modal scaling errors, characterized by a r , are disastrous for certain applications of experimental modal model, such 
as substructure synthesis, structural modification and adjustment finite element model^" 3 l Wei [4] analyzed the 
sources of modal scaling errors which are summarized as follows : 

1. Local and global calibration errors. 

2. Digital signal processing and FFT leakage errors. 

3. Improper orientation of the force or response transducer at the driving point. 

4. Low signal to noise ratio in the driving point measurement data. 

Another possible cause of modal scaling errors is non-linearities of structures. A typical nonlinear example given by 
Lallement^ showed that for a beam with nonlinear stiffness, lower natural frequencies and their corresponding mode 
shapes of a non-linear beam are very close to those of the same beam without the prescribed non-linear 
characteristics. However, a r could vary from 0.4 to 1.3 . 

From the above discussions, it can be concluded that modal scaling errors are almost inevitable in the applications of 
the prevailing modal testing techniques on practical structures. It is intended in this paper to determine the unknown 
scaler a r by a supplementary testing technique, i.e., mass additive technique. Some of the work in correcting modal 
scaling errors using mass additive technique^ has been done in the past using modal modeling technique. In this 
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paper, three analytical methods: first order and second order perturbation methods, and local eigenvalue modification 
technique are presented to predict the unknown modal scaling factors. 


FIRST ORDER PERTURBATION APPROACH 


The equation of motion for the general case of nonproportional damping is : 
My(t) + Cy(t) + Ky(t) = / (f) 
and its associated eigenvalue problem can be written as : 


A r 


c m r V'r 


M 0 


A r V»r 


-k a 


'Pr 


0 M 1 A r ^ r 


or in the condensed form : 

A r [U]{y r } = [V]{y r } 

where : 

'Pr 


<M-{A r * r }’ { ^ }eC2W 

[£/] = [£^|, [U]eR 2N , [U] = [Uf 

l y ] = In J- mei?” , [V] = [Vf 


( 1 ) 

( 2 ) 

(3) 


Defining a perturbated system : 

A r [U * eU]{1> r } = [V + eV]{^, T } (4) 

Where e is the perturbation parameter, and the elements of matrices V and U are of the same order as those of 
matrices [U] and [V] respectively. 

The eigenvalue A r may be expanded as power series in e : 

A r = A r +«A r ( 1 ) +e 2 A r W+ ••• (5) 

In the case of small perturbations ( e « 1 ), the higher order terms of e can be neglected, therefore Equation 5 can be 
written as : 

A r = A r +eA r ^) +0(e?) (6) 

The expression of A r (1) is given in numerous lectures f 7,8 ' : 

A r (1) = {y r } T [V - \ r U]{y r } (7) 

Equation 7 is based on the assumption : 
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{M T [tf]ty,}-i. (8) 

Considering the modal scaling errors, Equation 7 can be rewritten as : 

Ar (I) = {yrVlV - X r U]{y r W (9) 

Where {y T } is the identified modal vector with modal scaling errors: 

{y r } T [U]{y r } = VcZ ( 10 ) 

If the structural perturbation consists of only the mass matrix, AM = eM , Equation 9 can be written as : 

A r (1) = -W+rFm+r K 2 (11) 


In the present problem, the mass added structure is considered as the perturbated system and the perturbated 
eigenvalues A r and both eigenvalues and eigenvectors of the original system have been identified. Using small added 
masses, the unknown scaler a r can be directly derived from the first order approximation of A r : 

AA r = A r -A r = eA^ 1 ) = -A^, } T AM{i> r }a? (12) 


Thus : 


a, = v/AA^-A^FAM^,}) 

SECOND ORDER PERTURBATION APPROACH 


(13) 


It is important to note that the first order perturbation approach is applicable only in the case of small 
perturbations. However, from a practical point of view, it is desirable to hold AA r from 5% to 10% of A,. Therefore, 
the second order perturbation approach is introduced to improve the accuracy of the estimated modal scaling factors : 

AA r = A r -A r = eA r (1) +e 2 A r (2 > +0(e 3 ) (14) 

The expression of A r ® is : 

2N 

Ar (2) = E [-A r 2 {V'yF[MKV’r}a r a y ] 2 /(A r -Ay) - 2\ r {^ r } T [MMr}K {1) cZ (15) 

Obviously Equation 14 is nonlinear in e r . Therefore an iterative process is applied to calculate s r . This iterative 
process consists of three steps : 

1) Determining initial estimation of a r (0) : 

«r (0) = VU -A r )/(-A?{M r AM{lU) (16) 

2) Calculating : 

2N 

^A r < 2 )= S [-A r 2 {V’yF[AMKV-r} ar (0) aA0)]2/( Ar . A .). 2 A r 2 [{V. r nAMKV'r}(« r (0) ) 2 ] 2 (17) 

3) Recalculating : 
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a r W = y/(\ T -A r -eV^/^V’rfAW,}) 
Steps 2 and 3 are repeated until the convergence condition : 


(18) 


| a r * fc ) - a,**' 1 * | <a 

is satisfied by the k 01 iteration, where a is a given positive quantity. In general, a r (t) converges to a stable value very 
quickly. 


LOCAL EIGENVALUE MODIFICATION 


In most practical cases, the second order perturbation approach gives a reasonable estimation of modal scaling 
factors a r . But if the structure under study possesses very closed eigenvalues ( |(A r -A r+ i)/A r | <1% ), the 
perturbation method is no longer valid. Zhang and Lallement^ proposed a method based on modal space 
representation which would be suitable in this delicate case. This method requires not only the eigenvalues of the 
mass-added structure but also the modified eigenvectors. In this section, a method which does not require the 
eigenvectors of the mass added structure is presented : a method based on the local eigenvalue modification 
technique. 

Recall^ 10 ] that if the structure is perturbated only at a single degree of freedom, the eigenvalues of the perturbated 
structure are of the roots of the following equation : 

1/m = E(AA?#)/(A r -A) (19) 

r=L 

Where : 

A r and A are respectively the eigenvalues of the initial structure and the perturbated structure. 
rf> kr is the k component of the r* modal vector of the initial structure, 
m is the mass attached at the degree of freedom k. 

Because of the modal scaling error. Equation 19 becomes : 

1/m = £(AA?#a?)/(A r -A) (20) 

r=l 

Since the eigenvalues of the perturbated structure X 3 - ,j = 1,2, • • • ,N, are known, substitution of these values into 
Equation 20 yields : 


'1 Ini 


(Ai Ai^i)/ (Ai - Ai ) . . 

(Ai A n^In)/ &n “ ) 


1/m 


(Atf Ait/>Jti)/ (A x - X N ) . . 

(A N A n^In)I (^n ’^n) 

&N 


or in the condensed form : 

{Vm} = [Qlc?} (22) 

and then the unknown modal scaling factors { a? } are determined by solving Equation 22 : 

{a 2 } = [QYWrn] ( 23 ) 
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There are several comments needed to be made on the proposed local eigenvalue modification technique: 

1. Equation 20 implies the normalization relations : 

{y,mKM = i/«? . {y r } T [v][y r } = a m 

where : 

[^1= n. m=£o> > [V] = [V] T 

2. The major numerical difficulty in dealing with closely spaced modes is due to the following facts : 

• The difference between two adjacent eigenvalues A r - Ay appears in the denominators as shown in 
Equation 17. 

• It is difficult to make the correspondence between eigenvalues of the original structure and those of the 
mass-added structure. 

From Equation 21, it is noted that (i) the term ( A r - Ay ) does not appear in this equation, and (ii) the 
interchange of two modified eigenvalues does not affect the solution {a 2 }. Therefore, this method is very 
effective in dealing with closely spaced modes. 

3. The disadvantage of this method is that it is only applicable in the case of a single point and a single degree of 
freedom modification. In practice, this implies that the dynamics effect of the added mass must be dominant by 
only one of the six degrees of freedom at the mass attachement point. This theoretical restriction can cause real 
difficulty in choosing the proper location of the added mass. 

To avoid this problem, one of the options is to replace the additive mass by a feedback device , i.e., a device 
which can provide a force proportional to the acceleration at the mass attachment point i (/,• = m x , ). 
Another alternative is to modify the structure by adding a single spring instead of a mass. In case of adding a 
spring, Equation 20 becomes : 

■1/k = S(^a 2 )/(A r -A) (24) 

r= 1 

where k is the spring constant. 


PRACTICAL CONSIDERATIONS: Effective Mass 

Rotational degrees of freedom information at the mass attachment point(s) is also important in the application 
of the mass additive technique when the moment of inertia of the added mass(es) is not negligible. Various methods 
have been proposed to measure or predict the rotational degree of freedom frequency response functions or modal 
coefficients at a point on the structure^ 11,12 ], but all methods have their technical limitations. 

In this section, the concept of effective mass is introduced to include the dynamics effect of the unmeasured rotational 
degree of freedom information. Define the effective mass associated with the r 01 mode m efr as : 

m tit = m ad + A m r ( 25 ) 

where m ^ is the added mass and A m T represents the additional mass due to the effect of the rotational degrees of 
freedom of mode r. 

• 2 

The estimation of Am is based on the assumption that the kinetic energy due to rotation ( J © ) can be replaced by 


197 



. 2 

an equivalent translational kinetic energy (Ami ) along the z axis. 

For a one dimensional problem, total kinetic energy due to the added mass at point k can be written as : 
E = rriadzl +/0fc 

where J is the moment of inertia of the added mass. 

The angle of rotation © fc as shown in Figure 1 can be approximately calculated : 

@ k ~tangentG k = (z a 

where l ab is the distance between points a and b . 



Figure 1. Estimation of Rotational DOF from Translational DOF 


Equating A mi* to JQ k , i.e., E = ( rriaj, + A m ) z k , the scaler Am is determined : 

A T7l ~ J& k /z k = J (Z a z b ) / (JabZ-k ) 

Denoting E r as the total kinetic energy of the added mass associated with the r A mode : 
E r = A?(m ( *, + Am r )y& 

where : 

A m r = Jilpar-'Pbrf/dab'Pkr) 2 

Equation 28 shows that the effective mass m efr = m ^ + A m T is different for each mode. 


CASE STUDY 

In order to demonstrate the effectiveness of the proposed methods, a case study was performed on a *'g ht| y 
damped steel T plate shown in Figure 2. The Polyreference time domain modal parameter estimation method 1 was 
used to identify the eigenvalues and the eigenvectors of the unmodified T plate (4 references and 24 measurement 
points). Eigenvalues of the modified structure with two additive masses,i.e, mass 1: 97 grams at point 3 and mass 2: 43 
grams at point 22, were extracted from the frequency response functions measured at four reference and mass 
attachment points. The first mode shapes of the original T plate are plotted in Figure 3. 
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Figure 2. Test Structure: A Steel T Plate 






Figure 3. First Four Mode Shapes 


The first four damped natural frequencies and modal damping ratios of the original and the modified T plate are 
listed in Table 1. The identified complex modal vectors of the original T plate were normalized according to Equation 
8. The modal scaling factors a, calculated by using the first and second order perturbation methods are listed in Table 
2 . 

In order to verify the calculated modal scaling factors listed in Table 2, a structural modification software (DYNOP 
developed by U.C. SDRL) is used to predict the eigenvalues of the modified T plate with a single mass added at point 
3. Modified eigenvalues were derived for both the original and rescaled (using the modal scaling factors listed in Table 
2) experimental modal model. The predicted eigenvalues and the measured data are listed in Table 3. 
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TABLE 1. Measured Complex Eigenvalues 


Measured Complex Eigenvalues 


Mode No. 

Original Structure 

Modified T Plate(2 masses) 

f(Hz) 

?(%) 

f(Hz) 

<:(%) 

i 

177.37 

.984 

171.07 

1.322 

2 

334.20 

.441 

322.50 

.566 

3 

411.10 

.487 

392.02 

.558 

4 

581.7 

.286 

550.94 

.503 


TABLE 2. Calculated Modal Scaling Factors 


Calculated Modal Scaling Factors, a. 


First Order Approach 

Second Order Approach 

Mode No. 





real 

imaginary 

real 

imaginary 


1 

0.95 

0.05 



2 

1.05 

0.02 



3 

0.98 

0.04 

1.01 


4 

0.76 

-0.002 

0.68 

0.002 


TABLE 3. Complex Eigenvalues of the Modified T Plate 


Eigenvalues of the Modified T Platefl mass) 


Measured 

Predicted Eigenvalues 

Predicted Eigenvalues 


Eigenvalues 

Using Original {rp T } 

Using Rescaled {ip r } 

Mode No. 

f(Hz) 

<(%) 

f(Hz) 

Err(%) 

?(%) 

Err(%) 

f(Hz) 

Err(%) 

?(%) 

Err(%) 

HgjMB 

172.8 

1.289 

173.1 


.894 

-30.6 

172.8 

0 

1.166 

-9.5 


331.8 

.521 

333.2 


.469 

-10.0 

333.4 

.48 

.468 

-10.0 

1 1 9 1 1 

392.5 

.594 

389.7 


.284 

-52.0 

392.0 

-.13 

.483 

-18.7 


556.5 

.529 



.437 

-17.3 

563.9 

1.34 

.305 

-42.3 


From the results of this case study, it is noticed that: 

1. Since all modal testings of the T plate were very carefully conducted, and the structure is perfectly linear, 
therefore, all but the last a, are relatively small and within ±10% of 1.0 as shown in Table 2. 

2. It was pointed out in Reference[14] that the analytical predictions of modified eigenvalues using modal 
modeling technique could be corrupted due to modal truncation error. It was also indicated that the last few 
predicted damped natural frequencies are usually greater than the true ones. 

From Table 3, the 4 th damped natural frequency derived from the original modal vectors is less than the 
measured data. On the contrary, the one from the rescaled modal vectors is greater. This indicates that the 
damped natural frequencies based on the rescaled modal vectors are more reliable. 
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CONCLUSIONS 


The three analytical methods proposed in this paper can be used to improve the modal scaling factors of the 
measured modal vectors from the structure with or without close-spaced modes. The main advantage of these 
methods is that it only requires complex eigenvalues of the mass-added structure as additional information, i.e., only a 
few additional measurements are needed to obtain the complex eigenvalues of the mass-added structure. 

The concept of the effective mass is introduced to include the rotational effects of the additive mass(es) due to the 
lack of test data at the rotational degrees of freedom. Case studies based on this concept can be found in Reference 
[15]- 

The proposed methods also have the following limitations: 

1. The perturbation technique is applicable only in the absence of closely spaced modes. 

2. The local eigenvalue modification technique is limited to a single degree of freedom modification. 
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